Transit Time and Charge Correlations of Single Photoelectron 
Events in R7081 Photomultiplier Tubes 

F. Kaether, C. Langbrandtner 

Max-Planck-Institut fur Kernphysik, 
Saupfercheckweg 1, D-691 17 Heidelberg Germany. 
E-mail: Florian . Kaether@mpi-hd . mpg . de 



Abstract 

During the calibration phase of the photomultiplier tubes (PMT) for the Double Chooz experiment 
the PMT response to light with single photoelectron (SPE) intensity was analysed. With our setup 
we were able to measure the combined transit time and charge response of the PMT and therefore we 
could deconstruct and analyse all physical effects having an influence on the PMT signal. Based on 
this analysis charge and time correlated probability density functions were developed to include the 
PMT response in a Monte Carlo simulation. 



1 Introduction 

Double Chooz [IJ |2] is a reactor antineutrino ex- 
periment for the search of neutrino oscillations 
and a non vanishing neutrino mixing angle #13 
with two identical detectors located at different 
distances from the nuclear power plant at Chooz, 
France. Each detector will be observed by 390 
PMTs to record light pulses produced by neutrino 
induced interactions inside a gadolinium loaded 
liquid scintillator [3]. The reaction u e -\-p — >• e + +n 
leads to a coincident signature with e + annihila- 
tion and a delayed n-capture by gadolinium or 
hydrogen with an energy release of 8 MeV or 2.2 
MeV, respectively. With a light yield of 6000 pho- 
tons per MeV, a PMT coverage of 13% and a PMT 
sensitivity of 25% an average signal of approxi- 
mately 0.5 photoelectrons per PMT and MeV is 
expected. Therefore a precise knowledge of the 
characteristics of single photoelectron (SPE) sig- 
nals is absolutely necessary for the understanding 
of the Double Chooz detector behaviour. Exten- 
sive characterisations of all Double Chooz PMTs 
are described in [4j [5] . In this article we focus on a 
more detailed measurement of the PMT response 
concerning the charge and transit time character- 
istics of SPE events and their physical aspects. 
Furthermore, the obtained results were used to 



build probability density functions (PDFs) for the 
Monte Carlo simulations of the Double Chooz de- 
tectors and the data analysis. 

2 Experimental setup 

The hemispherical PMT R7081 (Hamamatsu Pho- 
tonics K.K.) has a diameter of 10 inch and a 
bialkali cathode. During the measurements, the 
PMTs are placed inside a Faraday dark room for 
light and electromagnetic noise protection while 
the complete electronics and DAQ are located out- 
side to avoid crosstalk. In addition, the PMTs 
are surrounded by a cylindrical mu-metal |6J to 
shield the terrestrial magnetic field. The shield- 
ing ensures a high energy resolution with a peak 
to valley ratio of ~ 4.5. As light source we use an 
A.L.S. picosecond laser with a wavelength of 438 
nm which is in the region of the PMTs highest 
sensitivity. The laser light is diffused and led by 
optical fibres to the front of the PMTs to create 
a complete illumination of the photocathodes. To 
create SPE signals we use a very low light inten- 
sity whereby only about 10% of the triggered light 
pulses cause PMT responses. While the number 
n of created photoelectrons is Poisson-distributed 
p(n) = \i n jn\ exp(— /i), the probabilities p(n > 1) 
are very low if \i « 0.1. 
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Figure 1: Schematic overview of the experimental setup. 



An overview on the electronical setup is shown 
in Figure [T] The Trigger system controls the laser 
output, defines the starting time for the Time 
to Digital Converter (TDC) and creates the time 
window for current integration with the Charge to 
Digital Converter (QDC). More details about the 
used electronical modules can be found in |4J. 

The value of the high voltage was calibrated 
for a gain of 10 7 . While the PMTs are back 
terminated with 50^ an average SPE charge of 
0.5 x 10 7 e = 0.8 pC is expected. The PMT sig- 
nals are decoupled from the high voltage with the 
splitter box. After amplification the charge is de- 
termined and a discriminator forms the stop sig- 
nal for the single-hit TDC. The VME based data 
acquisition allows to read out TDC and QDC si- 
multaneously for each event. 




Time [ns] 

Figure 2: Transit time spectrum of R7081 with SPE light 
intensity. 




Figure 3: Charge spectrum of R7081 with SPE light in- 
tensity. 



3 Analysis of the Time and Charge Spec- 
trum 

In this section we describe the effects occuring 
during the signal generation in the photomulti- 
plier tubes resulting in the observed structures in 
the transit time and charge distributions. 

Figure [2] and [3] show the transit time and charge 
spectrum of the R7081 photomultiplier using laser 
light with SPE intensity. Since the absolute value 
of the transit time is biased by cable lengths, pro- 
ceeding electronics etc. we focus on relative times 
while the maximum of the distribution is arbitrar- 
ily defined as t = 0. Regarding the relative tran- 
sit times in Figure [2j three types of events can be 
classified: 



• Main peak pulses: the major part of all PMT 
pulses (97%) are located in the time interval 
— 10 ns < t < 10 ns including mostly regularly 
multiplied electrons. 

• Pre-pulses: roughly 30 ns prior to the main 
peak so called pre-pulses occur, which arise 
from photoelectrons directly produced at the 
first dynode. 

• Late-pulses appear up to 70 ns after the main 
peak. They consist of photoelectrons which 
are (in-)elastically back scattered from the first 
dynode and produce secondary electrons after 
being reaccelerated to the dynode structure. 

Since the dark noise rate of the PMTs (caused by 
thermal emission of electrons from the cathode) 



Transit Time and Charge Correlations of Single Photoelectron Events in R7081 PMTs 



3 




Time [ns] 

Figure 4: Charge versus transit time distribution at SPE 
intensity. 
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Figure 5: Transit time of the Main peak fitted with two 
Gaussian curves and one expontential distribution convo- 
luted with a Gaussian with the same a as the central Gaus- 
sian. 

is in the range of a few thousand events per sec- 
ond, there is an observable accidental background 
in Figure [2] in form of a flat distribution with ap- 
proximately 5 events per bin. 

The charge of all events versus the transit time 
is shown in Figure [4j In the following sections 
we will investigate the different distributions and 
their physical origin in more detail. 

3.1 Main peak pulses 

Figure [5] zooms into the time distribution of the 
main peak in the time interval of —8 ns < t < 
10 ns. The observed structures can be described 
by the sum of two Gaussians (early and regular 
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Figure 6: Charge distribution of the main peak pulses. 
Black: Total charge spectrum. Blue: Regular pulses in 
a range of —2.5 ns < t < 2.5 ns fitted by two Gaussian 
functions. Green: Delayed pulses in the range of 2.5 ns < 
t < 8 ns fitted with a linear and a Gaussian functions. 
Red: Early pulses in the range of —6 ns < t < —2.5 ns 
fitted with two Gaussian functions. 



pulses, respectively) plus a convolution of an ex- 
ponential distribution with a Gaussian (delayed 
pulses) . In Figure [6] the charge distributions of 
the three different types of pulses are given, while 
their physical reasons are discussed in the follow- 
ing. 

Early pulses: The Gaussian distribution in the 
rising part of the main peak (see Figure [5| with 
a mean of Ie — —3.4 ns and ge = 1.0 ns includes 
the so called early pulses. Their origin can be 
explained by an elastic forward scattering of pho- 
toelectrons at the first dynode. These photoelec- 
trons keep their energy and directly hit the second 
dynode a few nanoseconds earlier than usual |8J. 
Because of the underlying effect one would expect 
the charge of these pulses to be reduced by the 
amplification factor of the second dynode which 
is estimated to be around 5. From the fit of these 
pulses with two Gaussians (see Figure [6| we get 
for the smaller one a mean charge of 0.15±0.02 PE 
which is in quite good agreement with an expected 
value of 0.18 PE. The second Gaussian around 
0.7 PE might be explained by inelastic forward 
scattered photoelectrons at the first dynode. 

Regular pulses: The central Gaussian with the 
mean at tn = — 0.2 ns and the so called transit 
time spread of cfr = 0.9 ns consists of the pre- 
dominant regular accelerated electrons. The tran- 
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sit time spread <jr is caused by different trajecto- 
ries of the photoelectrons through the accelerating 
electrical field mainly depending on their start- 
ing positions and velocities at the cathode. In 
contrary, during the multiplication process in the 
dynode system, a time-dispersion results in pulse 
broadening rather than in transit time spread be- 
cause of the large number of secondary electrons 
[9]. In Figure |6] the charge of the regular electrons 
is shown and described by a Gaussian around 
1.03 PE with a sigma of 0.3 PE. In addition, there 
is a small fraction of not optimal accelerated elec- 
trons falling into the applied time cut described 
by a second Gaussian at 0.3 PE and a sigma of 
0.14 PE. The relative probability ratio of both 
Gaussians is 20:1. 

Delayed pulses: The time distribution of the 
pulses in the falling part of the main peak is well 
described by a convolution of an exponential dis- 
tribution and a Gaussian (see Figure |5|. For the 
Gaussian a we expect the same underlying physics 
as for the transit time spread, therefore we choose 
o~d = cfr. The delay of these events can be ex- 
plained by photoelectrons produced at the outer 
side of the cathode running through a suboptimal 
acceleration due to field inhomogeneities between 
outer cathode and the first dynode. The latter will 
cause not only a delay but also a lower multiplica- 
tion at the first dynode, which can be seen in the 
charge distribution in Figure [6j where the charge 
fit is composed of a gaussian and an exponential 
(green). 

3.2 Pre- pulses 

Pre-pulses occur if the photon passes through the 
cathode and directly hits the first dynode. The 
difference of the transit time with respect to reg- 
ular pulses of roughly 30 ns is caused by the dif- 
ference of the electron time of flight in the ac- 
celerating field between cathode and first dynode 
compared to the speed of light. The charge of 
the pre-pulses is reduced by the amplification fac- 
tor of the first dynode which is expected to be 
around 15 for the R7081 series. Since the discrim- 
inator threshold for the PMT signal was 0.1 PE 
in our setup, we could only measure the higher 
part of the pre-pulse charge spectrum (Figure [7]). 
For this part we determined a pre-pulse probabil- 
ity of 0.1 %, but this number also depends on the 
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Figure 7: Charge distribution of the pre-pulses (blue) and 
whole spectrum (black). 

position of the light source and the angular distri- 
bution of the illumination. The events with charge 
values around 1 PE in figure [7] are accidental dark 
noise signals. 

3.3 Late-pulses 

Immediately after the main peak and up to 70 ns 
so called late-pulses occur. Since the nomencla- 
ture found in literature is inconsistent, we want 
to clearly distinguish late-pulses from after-pulses. 
The latter are created by the ionisation of resid- 
ual gas in the photomultiplier by photoelectrons. 
After-pulses have typical arrival times of some 
/is after the ionisation process and are there- 
fore always correlated to a previous PMT signal. 
Since the TDC measures only the next signal af- 
ter trigger-start, after-pulses do not influence the 
topic of this paper and are not further discussed. 
A description of the after-pulse behaviour of the 
PMT R7081 can be found in [10]. 

Late-pulses occur if the photoelectron is back 
scattered from the first dynode, then reacceler- 
ated by the electrical field and finally produces 
secondary electrons hitting the first dynode for the 
second time. This process of back scattering can 
happen elastic or inelastic producing heat at the 
dynode |7J. This explanation is consistent with 
the observation that the arrival time of late-pulses 
ends at around 70 ns after the regular pulses which 
is approximately twice the time of flight of the 
photoelectrons between the cathode and the first 
dynode (= pre-pulse time). It is slightly longer 
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Figure 8: Charge distribution of the late-pulses. Black: 
total charge spectrum. Blue: charge spectrum of 
mostly elastic scattered late-pulses with 60ns<t<75ns. 
Red: charge spectrum of inelastic late-pulses with 
15 ns< t <60 ns, both fitted by a linear and a Gaussian 
function. 



since the path length of the back scattered and 
back accelerated electron is not a straight line but 
rather a curve in most cases. 



To distinguish between inelastically and elas- 
tically scattered photoelectrons we separate the 
late-pulses by their arrival times before and after 
60 ns, respectively. This separation is motivated 
by the expectation of shorter path lengths of in- 
elastically scattered photoelectrons. The charge 
distibutions of both samples are shown in Figure 
[8j The separation of (in-) elastic events does not 
work perfectly well, however the fraction of events 
with lower energies is higher at earlier times, as 
expected. The charge distributions are described 
by a linear and a Gaussian part, while the origin 
of the linear distribution is expected to be caused 
by inelastically scattered electrons. The Gaussian 
charge distribution is expected to be generated by 
elastically scattered photoelectrons. 

In our setup we measured the probability for 
late-pulses after t>10ns to be 3%. In the de- 
scription of the Monte Carlo simulation given in 
the next section we will see that an elastic to in- 
elastic late-pulse ratio of 1:2 fits the data best. 



Pulse 


Time 


Charge 


Prob. [%] 


Pre-pulses 


G 


G 


0.1 


Early-pulses 


G 


2G 


1 


Middle-pulses 


G 


2G 


76.9 


Falling edge pulses 


exp 


exp 


19 


Inelastic late-pulses 


exp 


exp 


2 


Elastic late-pulses 


G 


G 


1 



Table 1: PDFs for the different pulses. G stands for a 
Gaussian, exp for an exponential probability distribution. 



4 Simulation of the PMT response 

4.1 Probability density functions 

In the Double Chooz experiment the PMT charge 
and time response influences not only the energy 
reconstruction of detector events (e.g. from neu- 
trino interactions), but also higher order analy- 
sis tools as vertex reconstruction and pulseshape 
analysis. Therefore it is crucial to understand the 
transit time and charge behaviour of the PMTs 
and to include them into the detector simulation 
in an appropriate way. In the previous section it 
was shown that effects inside the PMT having an 
influence on charge and time response can be well 
described by six different types of pulses, namely 
one type for pre-pulses, three types for main peak 
pulses and two types for late-pulses. The mea- 
surements described in the previous sections are 
the basis for probability density functions (PDFs) 
which are derived from the data fits for all these 
types of pulses and are listed in Table 1. Their 
composition in transit time and charge is given in 
Figure [9j 

4.2 Monte Carlo Simulation result 

The steps of the Monte Carlo Simulation are as 
follows: first it is decided by a Poisson distribution 
how many photoelectrons are created at the cath- 
ode. The next step decides which type of pulse 
from table 1 occurs and according to the PDFs 
for this pulse type the charge and transit time is 
chosen. This step is done for each photoelectron 
separately. The last step calculates the sum of all 
occurring photoelectrons which finally gives the 
total charge of the event while the transit time is 
calculated by the shortest arrival time of all pulses 
in the event. 

The probabilities of different pulses were cho- 
sen to be in best agreement with the measured 
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Figure 11: Simulated (red) and measured (black) time (left) and charge (right) distribution for a multi PE light intensity 
with a Poissonian mean of \i = 7.7. 
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Figure 12: Simulated charge versus transit time distribu- 
tion. 



time and charge spectra. The results can be seen 
The charge versus transit time for 



in Figure 10 



all simulated events is shown in Figure [12] and 
can be compared with the measured distribution 
in Figure [4] 

Also for higher light intensities the simulation 
yields good results compared to the measured dis- 
tributions which is illustrated in Figure [TT] for an 
incident mean number of photoelectrons at the 
photocathode of /i = 7.7. In the real transit time 



spectrum of figure [TT] the pre-pulse and main peak 
distribution is wider than in the simulated spec- 
trum, as for this measurement not a laser source 
but a LED with a larger timing jitter was used. 
Furthermore, a forward shifting of the main peak 
of 2 ns can be explained by the fact that larger 
pulses arrive earlier at TDC threshold. This effect 
is not caused by the PMT but by applied electron- 
ics (the use of a leading edge discriminator) and 
is therefore not included in the simulation. 

5 Summary 

The single photoelectron response of PMT R7081 
was analysed regarding charge and transit time 
correlations and explained by the underlying 
physical effects. The results of the measurements 
were used to build probability density functions 
for a usage in Monte Carlo simulations which were 
finally implemented into the Double Chooz detec- 
tor simulation software. 
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